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Scattering of coherent light as it propagates parallel to a shock wave, formed in front of a bluff 
cylindrical body placed in a supersonic stream, is studied experimentally and numerically. Two 
incident optical fields are considered. First, a large diameter collimated beam is allowed to pass 
through the shock containing flow. The light intensity distribution in the resultant shadowgraph 
image, measured by a low light CCD camera, shows well-defined fringes upstream and downstream 
of the shadow cast by the shock. In the second situation, a narrow laser beam is brought to a grazing 
incidence on the shock and the scattered light, which appears as a diverging sheet from the point of 
interaction, is visualized and measured on a screen placed normal to the laser path. Experiments are 
conducted on shocks formed at various free-stream Mach numbers, A/, and total pressures, P 0 . It is 
found that the widths of the shock shadows in a shadowgraph image become independent of M and 
P 0 when plotted against the jump in the refractive index, An, created across the shock. The total 
scattered light measured from the narrow laser beam and shock interaction also follows the same 
trend. In the numerical part of the study, the shock is assumed to be a “phase object,” which 
introduces phase difference between the upstream and downstream propagating parts of the light 
disturbances. For a given shape and An of the bow shock the phase and amplitude modulations are 
first calculated by ray tracing. The wave front is then propagated to the screen using the Fresnel 
diffraction equation. The calculated intensity distribution, for both of the incident optical fields, 
shows good agreement with the experimental data. © 1995 American Institute of Physics. 



I. INTRODUCTION 

A significant change of gas density across a shock is 
accompanied by a change in the refractive index. Upstream 
of the shock is a region of lower refractive index and down- 
stream is a region of higher refractive index. A large deflec- 
tion of light rays across this change in refractive index makes 
the shock easily visible in schlieren and shadowgraph im- 
ages. The phenomenon is usually explained with the laws of 
geometrical optics (Merzkirch, 1 p. 115). However, a shock is 
also expected to cause light diffraction, as the extremely 
small thickness effectively makes it a sharp interface be- 
tween the optically different upstream and downstream re- 
gions. The objective of the present work is to explore the 
light diffraction effects caused by a shock wave. This paper 
presents extensive experimental evidence and a numerical 
calculation based on the Huygens-Fresnel principle 2 to ex- 
plain this phenomenon. 

Two different optical conditions, both of which involve 
significant diffraction effects as laser light propagates paral- 
lel to a shock surface, are considered. The first one is the 
well-known shadowgraph situation, where a large diameter 
collimated laser beam is allowed to pass through the shock 
containing region and the shadowgraph image formed on a 
screen is studied. 

The second is somewhat unknown and was first reported 
by Panda. 3,4 It involves a narrow laser beam brought to a 
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grazing incidence on a shock surface. In this situation a thin 
diverging light sheet is found to appear from the point of 
interaction. The visible divergence angle of the light sheet is 
very large, up to ±10° depending upon the shock strength, 
and the light appears both upstream and downstream of the 
shock. When a screen is placed normal to the laser path, the 
cross-section of the sheet appears as a long, bright streak. 
The scattered light is found to disappear at all other inci- 
dence angles when the laser beam pierces the shock. The 
optical phenomenon is also believed to appear primarily due 
to diffraction of laser light by shock waves; however, suffi- 
cient evidence to determine the exact nature of the phenom- 
enon was absent. It has been demonstrated that the phenom- 
enon can be used as the basis of a novel shock detection 
technique 3,5 which is convenient in detecting unsteady shock 
motion. A major motivation of this work is to determine the 
effects of various fluid dynamic and optical parameters on 
the scattered light, as well as to establish the physics behind 
the phenomenon. 

Evidence of diffraction fringes in the shadow image of a 
shock can be found in many earlier experiments. 6,7 Pfeifer 
et a/. 8 have shown that for a plane shock, the zeroth diffrac- 
tion order is a minimum, which appears in the shadowgraph 
to form the dark band. Merzkirch 1 also discusses that, when 
the incident light waves propagate parallel to the surface of 
the shock the major contribution to the formation of a 
shadow is light diffraction (see discussion on p. 133). How- 
ever, all of these discussions are very limited in scope and 
frequently the diffraction effects are overlooked. Wide angle 
scattering of a narrow beam of light by conducting cylinders 
has been reported earlier in the literature (see Langois 
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et al 9 ). However, the present problem involving shock 
waves is very different. Optical properties of the scattered 
light such as the intensity distribution, fringe spacing, and 
the effects of various fluid dynamic parameters were not 
studied for the diffraction effects caused by a shock wave. 

In the experimental part of this study a steady bow 
shock, formed ahead of a bluff cylinder placed in a super- 
sonic free stream, is used to study such effects. Light inten- 
sity distributions for both the wide angle scattering from a 
narrow beam and diffraction fringes in shadowgraph images 
are measured for various free-stream Mach numbers and to- 
tal pressures. In the numerical part of this study, the phenom- 
ena are modeled as due to the phase difference between the 
light waves traveling upstream and downstream of the shock. 
The two-dimensional (2-D) analysis assumes that the inci- 
dence disturbance field is made of plane and parallel waves 
to model the shadowgraph image and a Gaussian strip to 
model the narrow laser beam and shock wave interaction. 
The diffraction patterns are calculated and presented after 
describing the experimental results. Finally, the similarities 
and differences between the two situations and various other 
parameters that may affect the diffraction phenomena are 
discussed in the conclusion. 

II. EXPERIMENTAL FACILITY 

The experiments were conducted in a continuous, 127 
mm XI 27 mm supersonic wind tunnel at the NASA Lewis 
Research Center. Details of the tunnel can be found 
elsewhere. 10 In brief, the tunnel is supplied by compressed 
air which is discharged into a low- vacuum altitude exhaust 
system. The free-stream Mach number can be varied from 
1.5 to 3.0 by changing the nozzle contour. The free-stream 
total pressure, P 0 , also can be varied independently from 134 
to 375 KPa. The total temperature of the supplied air is the 
same as that of the ambient air (294 °K). For the present 
experiment the test section was equipped with two scratch 
free plexiglass windows, each 25.4 mm thick. The cylindri- 
cal bluff body used to produce a bow shock was mounted 
edge on with its axis along the flow direction. It was 12.7 
mm in diameter and was mounted on a sting at the center of 
the tunnel. 

The experiments were conducted at three different free- 
stream Mach numbers: 1.5, 2.(3, and 2.5, and five different 
free-stream total pressures: 375, 306, 237, 168, and 134 KPa. 
For the highest Mach number of 2.5 the tunnel could not be 
started at the lowest total pressure of 134 KPa, and for the 
lowest Mach number of 1 .5 operations were limited to a 
maximum P {) of 237 KPa by the maximum allowable mass 
flux of air supplied to the tunnel. 

A 4 mW, random polarization. He-Ne Laser (632.8 nm 
wavelength) was used as the light source for all experiments. 
A polarizer, placed in front of the laser was used to change 
incident beam polarization, and a combination of the polar- 
izer and a neutral density filter was used to cut down the 
laser power when needed. 

The scattered light patterns were visualized on a semi- 
transparent screen (similar to a ground glass) placed normal 
to the laser path. The images on the screen were photo- 
graphed using a cooled CCD camera which provided quan- 



titative information of the light intensity distribution on the 
screen. The 12-bit, 576X384 pixel camera has a separate 
refrigeration unit to cool the camera head to -45 °C; thereby 
dropping the dark current significantly. This low noise fea- 
ture makes the camera useful for measuring low levels of 
light. The camera has its own controller which, in turn, is 
interfaced with a 486DX Personal Computer (PC). Images 
obtained by the CCD camera were stored in the PC and are 
post-processed to obtain quantitative information. Various 
Nikon lenses were used on the camera head to expose a wide 
range of area on the screen. 

A laser shadowgraph system was used to obtain a large 
diameter collimated beam (Fig. 1). The narrow beam from 
the He-Ne laser was allowed to pass through a beam ex- 
pander (made up of a microscope objective, a 25 gm spatial 
filter and a 300 mm focal length lens) and was allowed to 
pass through the test section. The shadowgraph image 
formed on the semitransparent screen placed on the other 
side of the test section was photographed either by the CCD 
camera or by a 35 mm camera. 

III. EXPERIMENTAL RESULTS 

Figure 2 shows shadowgraph images of the bow shock 
formed in front of the bluff body for two different free- 
stream Mach numbers. For all figures presented in this paper, 
flow is from left to right. The bow shock is found to become 
more curved and to move closer to the bluff cylinder as the 
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FIG. 2. Shadowgraph image of the bow shock formed ahead of a bluff 
cylindrical body; (a) M = 1.5 and (b) 2.0. Flow is from left to right. 
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free-stream Mach number is increased. However, for a fixed 
A/, as P {) is changed the shock shape and its location are not 
affected. 

A steady shock is essential to achieve the experimental 
goal. To check the shock steadiness, video recordings of the 
shadowgraph images were obtained by a separate light inten- 
sified CCD camera that can be gated to a fast shutter speed of 
20 ns. A frame by frame analysis showed that the shock was 
extremely stationary up to A/^2.5, above which it mean- 
dered at a very low frequency, perhaps due to the proximity 
to the bluff body. 

A. Diffraction effects in shadowgraph image of shock 

The diffraction effects are not visible in the 35 mm pho- 
tographs of Fig. 2. To highlight the detail of the shock 
shadow a small region around the tip of the bow shock along 
the centerline of the cylindrical body was photographed us- 
ing the low light CCD camera. Two different images were 
taken: one without the any airflow, “no-flow" image, and a 
second one with the flow on. Subsequently, the no-flow im- 
age was subtracted from the “flow" image to determine the 
changes caused by the shock. The processed data clearly 
show the dark shadow cast by the shock and clear fringes 
upstream and downstream of the shock. Figures 3(a) and 
3(b) show two such processed images of the bow shock in 
M- 1.5 free-stream at two different free-stream total pres- 
sures. The screen, from which the shadowgraph images were 
photographed, was placed 685 mm away from the shock lo- 
cation. The dark band representing the shadow of the shock 
appears very thick due to a magnification of —8 from Fig. 2. 

A special problem encountered with the laser shadow- 
graph system is the speckle noise produced by small dust 
particles on the wind tunnel windows. This random optical 
noise is easily visible in Figs. 3(a) and 3(b). To cancel out the 
background noise the light intensity data was averaged over 
a transverse band of 30 pixels. The intensity values presented 
in Fig. 3(c) are also normalized by the average intensity of 
the no-flow image. The resultant intensity distribution shows 
that the first fringe downstream of the central dark band is 
always the brightest. Both downstream and upstream fringes 
are indicative of light diffraction, whereas some of the extra 
light in the downstream fringes can also be interpreted from 
geometrical optics due to refractive bending of light. 

It is found that the width of the central dark band and the 
intensity modulation in the fringes depend on the free-stream 
total pressure P 0 and Mach number A/. For a fixed M, the 
light patterns become weaker as P {) is decreased. This is 
clearly visible from the two different P {) cases shown in Fig. 
3. It is possible to simplify this dependency by analyzing the 
fluid dynamic and the optical properties of the shock. Table I 
presents such properties calculated for the tip region of the 
bow shock. Since the shock is normal to the flow direction at 
this location, flow properties downstream of the shock are 
obtained by applying the normal shock relations. The refrac- 
tive indices (n) are calculated assuming ideal gas behavior, 
and by applying the Gladstone-Dale relationship: 1 n~ 1 
= Cp, where p is the local density of air and C is the 
Gladstone-Dale constant (C=2.25Xl(T 4 m 3 /Kg used for 
the laser wavelength). Data in Table I show that as P 0 is 



Distance from cylinder (mm) 


FIG. 3. (a) and (b) Enlarged view of the tip region of the bow shock ob- 
tained with a CCD camera; M ~ 1.5, P 0 =237, and 134 KPa, respectively, (c) 
Light intensity distribution along centerline, — P 0 =237 KPa and 
-_--/> (| =]34 KPa. The latter is shifted by one major division. 


decreased keeping M constant, the air density, both upstream 
and downstream of the shock, decreases. The refractive indi- 
ces also follow the same trend. The last column of Table I 
shows that the net difference of the index of refraction, A/a, 
between the upstream and downstream regions decreases 
with a decrease in P {) . The dependence on M for a fixed P () , 
however, is not straightforward as the density change across 
a shock is a nonlinear function of M . Nonetheless, it is found 
that the jump in the refractive index across a shock is the 
single most important parameter that controls the shadow 
formation in the shadowgraph images. This is shown in Fig. 
4 where the width of the central dark band, representing the 
shadow of the shock, is plotted against A n. The data points 
are obtained from shadowgraph images on a screen placed at 
a fixed distance of 216 mm away from the shock and for 
different M and P 0 conditions. The normalized light inten- 
sity distributions, similar to Fig. 3(c), were obtained from 
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TABLE I. Shock properties relevant to the light scattering phenomenon. 


Mach No. M 

Freestream 
total 
Pressure 
P 0 , KPa 

Air density 
Kg/m 3 


Refractive index 


Upstream Downstream 

Upstream 

Downstream 

Up, n u 

Down, n d 

Diff., An 


237 

1.115 

2.076 

1.000 252 

1.000 468 

0.000 217 

1.5 0.70 

168 

0.791 

1.473 

1.000 178 

1.000 332 

0.000 154 


134 

0.629 

1.171 

1.000 142 

1.000 264 

0.000 122 


375 

1.027 

2.738 

1.000 232 

1.000 618 

0.000 386 


306 

0.838 

2.235 

1.000 189 

1.000 504 

0.000 315 

2.0 0.58 

237 

0.649 

1.732 

1.000 147 

1.000 391 

0.000 244 


168 

0.461 

1.228 

1.000 104 

1.000 277 

0.000 173 


134 

0.366 

0.976 

1.000 083 

1.000 220 

0.000 138 


375 

0.588 

1.959 

1.000133 

1.000 442 

0.000 309 


341 

0.534 

1.779 

1.000 12 

1.000 401 

0.000 281 


306 

0.48 

1.6 

1.000 108 

1.000 361 

0.000 253 

2.5 0.51 

272 

0.426 

1.419 

1.000 096 

1.000 320 

0.000 224 


237 

0.372 

1.239 

1.000 084 

1.000 280 

0.000 196 


168 

0.264 

0.879 

1.000 059 

1.000 198 

0.000 139 


each image and the thickness of the shadow is measured as 
the width between the unit intensity points of the trough 
representing the shock shadow. The refractive index jump, 
A n y is obtained from Table I. Figure 4 shows that the shadow 
width increases with An and is independent of M and P 0 . 
The implication of this result is that the only way any change 
in the fluid dynamic parameters affects the optical scattering 
is by changing the optical property of the medium, that is the 
net jump in the refractive index across a shock. Somewhat 
similar behavior is also observed when a narrow laser beam 
is used to study the shock laser interaction. This is discussed 
next. 

B. Narrow laser beam and shock interaction 

Figure 5 shows a schematic of the arrangement to mea- 
sure the scattered light produced from such an interaction. 
When the narrow laser beam is brought to a grazing inci- 
dence on the bow shock, the scattered light appears as a 
diverging sheet from the main laser beam. For all data pre- 



Refractlve Index Jump (An) 


FIG. 4. Width of the central dark band in shadowgraph image of shock 
versus the refractive index difference created across the shock; •, M = 1.5; 
▲, M ==2.0; ■, A/ = 2.5. 


sented in this paper, the laser beam is made to touch the tip 
region of the ‘bow shock. The light patterns, as photographed 
from a screen placed normal to the laser path and 480 mm 
away from the shock, are shown in Fig. 6. The long streaks, 
extending far upstream and downstream from the shock, are 
aligned with the flow direction. If the laser beam is moved to 
touch any other part of the bow shock, the streaks are found 
to become inclined along the local normal to the shock sur- 
face. Various other features including the directional depen- 
dence of the scattered light has been discussed earlier. 3,4 The 
scale at the bottom of the figure shows the angular spread (6) 
of the scattered light. The angular spread angle is defined as 

f^tan ^ 1 

/ Distance on screen from center of the laser spot 

\ shock to screen distance 

(1) 

and is also shown in the schematic of Fig. 5. The wide angle 
light streaks are definite evidence of diffraction effect. Fig- 
ures 6(b)- 6(d) correspond to a decreasing P Q while the free- 
stream Mach number was kept constant. Clearly, the phe- 
nomenon weakens progressively as P 0 is decreased. Note 
that the laser beam, whose cross section is seen as a bright 
central spot, splits into two parts after touching the shock. 


Binrr 

Cylinder 



Coofttd CCD 

Camera 



FIG. 5. Schematic of the arrangement to visualize and measure the wide 
angle light scattering phenomenon. 
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FIG. 6. Cross-sectional views of the laser beam showing the wide angle 
light scattering; (a) no-flow; (b)-(d) interaction with shocks formed in 
M =2.0 flow and P o =306, 237, and 168 KPa, respectively. 


Quantitative measurement of the light intensity distribu- 
tion in the long streaks was not straightforward, as the inten- 
sity varied by many orders of magnitude from the center of 
the laser beam to the visible end of the streaks. To satisfy the 
conflicting requirements of avoiding saturation from the in- 
tense main laser beam and having sufficient exposure to cap- 
ture the weak tail, data were obtained in two steps. At first, 
the laser power is reduced considerably by placing a neutral 
density filter in the path of the beam, and the paraxial region 
around the main beam is photographed. The images show 
only the change in the main laser beam; little information 
about the scattered light is present. Figures 7(a) and 7(b) 
show two such photographs of the main laser beam with and 
without the presence of the scattering phenomenon. The cen- 
tral dark region between the split portions of the beam in Fig. 
7(b) corresponds to the shadow formed by the shock. The 
light intensity distribution measured from the pixels along 
the image centerline is plotted in Fig. 7(c). The intensity data 
are normalized by the maximum value in the undisturbed 
laser beam. Clear splitting of the beam and spreading of light 
energy primarily in the downstream direction are visible in 
this figure. A comparison of the scales showing the spread 
angle in Figs. 6 and 7 shows that the latter corresponds to a 
small region close to the main beam. 

A second set of images was taken to measure the inten- 
sity distribution in the wide angle streaks (Fig. 8). For these 
measurements the incident laser beam had to be brought 
back to its full intensity and a 12 mm wide beam stop was 
placed just in front of the semitransparent screen to block the 
main beam. The latter was necessary to avoid blooming of 
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FIG. 7. Enlarged view of the main laser beam; (a) no-flow; (b) interaction 
with a shock in 2 and P {) ~ 237 KPa flow; (c) light intensity distribution 
along centerline; , no-flow; — M=2.0, P 0 = 237 KPa. 

the CCD elements from the intense light of the full power 
beam. The shadow cast by the beam stop is marked in this 
figure. The diffraction effects from the beam stop itself are 
negligible as it is placed just in front of the screen. Once 
again, a comparison between the abscissa of Figs. 7 and 8 
shows that the entire paraxial region covered in the former 
falls into the shadow of the beam stop. 

The intensity distribution in the long light streaks for 
various upstream total pressure conditions is shown in Fig. 8. 
The free-stream Mach number is kept constant at 2.0 for all 
of these measurements. All data are normalized by the maxi- 
mum intensity value measured in the M=2.0, F 0 =375 KPa 
condition. The intensity of the scattered light is found to 
diminish with a drop in the free-stream total pressure, P {) . 
Recall that P () determines the jump in the refractive index 
across the shock. The scenario is consistent with the earlier 
description of shadowgraph images. 

For an application of the light scattering phenomenon as 
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FIG. 8. Intensity distribution in the scattered light from interaction with 
shocks formed in M =2.0 flow. Shock to screen distance=292 mm. 


a detection technique for shock, it is useful to obtain a rela- 
tionship between the net scattered light and the shock 
strength. The latter is measured as the jump in the refractive 
index, An, across a shock. The net scattered light, / r , is 
expected to be proportional to the area covered under each 
curve of Fig. 8. A plot of I T versus An is shown in Fig. 9. All 
data were obtained from measurements with a fixed laser 
power and a fixed shock to screen distance. 

The close-to-linear dependence in the semi log plot of 
Fig. 9 shows that I T increases almost exponentially with An. 
Unlike the shadowgraph situation, I T is also dependent on 
the free-stream Mach number. However, the calculated val- 
ues of l T are not very accurate, as the scattered light from the 
paraxial region of beam propagation was blocked by the 
beam stop, and therefore, is excluded from the calculation. 

The effect of the polarization of light on the wide-angle 
scattering phenomenon was studied by placing a polarization 
rotator on the path of the incident laser beam. However, 
within the measurement accuracy, the phenomenon is found 
to be independent of the light polarization. 
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FIG. 9. Net scattered light in narrow laser beam and shock interaction 
versus the refractive index jump across the shock; •, Af = 1.5; A, M =2.0; 
■. M = 2.5. 


IV. ANALYSIS 

According to the laws of geometrical optics, light rays 
are expected to bend only along the gradient of refractive 
index, which is along the downstream direction for a normal 
shock. However, clear violations are observed in the pres- 
ence of the upstream fringes in the shadowgraph image and 
in the upstream part of the scattered light in narrow laser 
beam and shock interaction. The deflection angle, y, of a 
laser beam while crossing a shock has been calculated by 
Kriksunov and Pliev 11 following the paraxial ray equation, 
tan y=/o(</rt/d£)dz, (where z is the beam propagation di- 
rection and £ is the flow direction) and is found to be a 
fraction of a degree. The visible spread angle of the long 
streaks, on the other hand, can be as large as ± 10°. 

Shock waves are also known to cause light reflection. 
However, the reflectivity of the shock fronts are very small. 
Measurements by Homig 12 showed that the optical reflectiv- 
ity of a strong shock is in the order of 10~ 6 , while the present 
data show that the net scattered light from a narrow laser 
beam and shock interaction can be greater than 1% of the 
incident laser power. Therefore, contribution from the re- 
flected light to the net scattered light is expected to be mini- 
mal. 

The light scattering phenomenon in the two situations 
considered here are indicative of a strong diffraction effect. 
Noticeably, both phenomena appear when light disturbances 
travel parallel to the surface formed by a shock. In this con- 
dition, a part of the beam propagates upstream and the rest 
downstream of the shock. The difference in the optical path 
length between the two parts of the beam produces phase 
variation. Therefore, a shock wave can be thought of as a 
“phase object” 2 that does little change to the magnitude of 
the incident light disturbance, but considerably distorts the 
phase. The resulting diffraction pattern can be calculated by 
an application of the Huygens -Fresnel principle. In an el- 
ementary way, Pfeifer et a/. 8 used the above principle to cal- 
culate the light intensity variation caused by the passage of a 
normal shock in a shock tunnel. It may be worthwhile to 
mention that any compressible flow, in general, is a phase 
configuration due to the variation in the refractive index with 
little absorption taking place, and a shock wave is the ex- 
treme case. 

The present analysis is performed for the two incident 
light disturbance fields for which experiments were per- 
formed. At first, the ray tracing method is used to calculate 
the phase and intensity variation in a bundle of light rays just 
after crossing the shock wave. Next, scalar diffraction theory 
is applied to calculate the irradiance distribution after the 
emerging rays propagate to the screen. 13 Figure 10 shows a 
schematic of the steps involved. The calculation is performed 
in 2D and for the tip region of the bow shock. The £ axis is 
along the flow direction and the z axis is along the light 
propagation direction. The shape of the bow shock is ob- 
tained from the shadowgraph image. A simple curve fit, 
€=a[l — cos (bz)], where a and b are adjustable constants, is 
found to be sufficient to analytically express the shock pro- 
files. The values of the constants for the M = 1.5 and 2.0 
cases are, respectively, a~ -1-68.62 and +70.5 mm, and 
b =0.01 91 and 0.0232 mm -1 . The refractive indices are as- 
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sumed to be constant, n u , everywhere upstream, and n d , 
everywhere downstream of the shock. 

The ray tracing calculations start with an initial intensity 
distribution 0 O (£) fr° m position z — — Zo > away from the 
shock. The deflection of each ray, while crossing the shock 
front, is determined by an application of Snell’s law. For 
example, at point 1 in Fig. 10, sin SJsin 8 d~ n J n u » where 
is the incidence angle and 8 d is the refracted angle. After 
each ray emerges from the shock it is traced back to a virtual 
line of origin, z — z v , to avoid ray crossing and also to obtain 
a standard reference line to start the diffraction calculation. 
The point of intersection of the traced back first ray, above 
£=0, to its original path is used to determine z v . The com- 
plex wave front at this virtual line of origin is determined as 
^i(4)=a(4)exp[i^(£)], where phase $£)=fcS(£), where S(£) 
is the computed optical path and k is the wave number 27 r/X; 
X is the wavelength of the laser light. The amplitude a(g) is 
determined using conservation of energy along ray bundles. 
Assuming the rays do not cross, the energy between 0 to £ 0 at 
the starting line should be the same as that between 0 to £ at 
the virtual line of origin. This provides the integral equation 

f ip 0 ip$ [ tpjipf d £, ( 2 ) 

Jo Jo 

where the superscript * represents the complex conjugate. 

In the next step, the calculated wavefront is propa- 
gated to the screen, z distance away from the shock, using 
the Fresnel diffraction equation, 

Ke ikz f* 2 

Ip(x)= tp,(i)e {ik ' 2z)ix - f) d£, (3) 

Z J —00 

where ip(x) is the optical disturbance on the screen and K , 
the “obliquity factor ” = — i/X. The complete procedure to nu- 
merically calculate the Fresnel integral is provided by 
Weaver 14 (see Chap. 8) and will not be repeated here. In 
brief, it involves casting the above equation in a convolution 
form and applying Fourier transform and an inverse trans- 
form. The calculations use Fast Fourier Transform routines 
from the Numerical Recipes 15 and are performed in the same 
PC used for data acquisition. 

Figure 11(b) shows the calculated light intensity distri- 
bution on a screen when the incident light disturbance 0 0 (£) 
is assumed to consist of plane, parallel waves. This situation 
corresponds to the shadowgraph image of the shock. For a 



Distance from shock (x, mm) 

FIG. 11. A comparison between (a) experimental and (b) calculated inten- 
sity distribution in the shadowgraph image of shock; M - 2.0, P 0 ~231 KPa, 
shock to screen distance=216 mm. 


direct comparison, experimental data obtained from an iden- 
tical condition is shown in Fig. 11(a). All major features, 
such as the central shadow and the fringe patterns upstream 
and downstream of the shock, are clearly captured in the 
numerical calculations. The width of the shock shadow in the 
experimental and the numerical data are nearly the same. 
However, the fringe spacings are slightly different and the 
peak intensity in the fringes are higher in the calculated data. 
The latter is due, possibly, to an error in the experimental 
data, which is contaminated by a high level of speckle noise. 
Moreover, as mentioned earlier, the experimental data were 
also averaged over 30 pixels. 

Note that the zero on the abscissa of Fig. 1 1 represents 
the position of the shock, positive numbers represent down- 
stream and negative numbers represent upstream distance. 
The numerical data in Fig. 11(b) demonstrate that the shock 
is located at the upstream edge of the thick shadow. This is 
consistent with the discussion of Merzkirch. 1 

It is interesting to note an important difference between 
the present analysis of light scattering by shock and the clas- 
sical Fresnel diffraction by a straight edge. 2 Both problems 
are similar in the presence of a sudden discontinuity; how- 
ever, there is an important difference between the two. The 
textbook treatment of diffraction by a straight edge assumes 
an opaque half-screen that completely absorbs the incoming 
light, while a shock wave does not cause any light absorp- 
tion. For the latter, transmitted light beyond the “straight 
edge” created by a shock produces the downstream fringes in 
the shadowgraph situation, and for the former, light intensity 


Phys. Fluids, Vol. 7, No. 9, September 1995 


J. Panda and G. Adamovsky 2277 




FIG. 12. A comparison between (a) experimental and (b) calculated inten- 
sity distribution in the narrow laser beam and shock interaction; M~ 2.0, 

P 0 =237 KPa, shock to screen distance = 737 mm; undisturbed beam; 

— after interaction with shock. 


approaches zero very quickly inside the geometrical shadow 
of the opaque screen. 

To calculate the light scattering from the narrow laser 
beam and shock interaction, the incident light intensity dis- 
tribution is changed to that of a 2 mm diameter Gaussian 
beam with a peak intensity of unity, ^(£)=exp(— £ 2 /4). The 
beam diameter corresponds to the estimated experimental 
value at the shock position. Since the calculations are per- 
formed in 2D, the incident field is, in effect, a Gaussian strip 
of light. The calculated diffraction pattern is shown in Fig. 
12 along with the experimental data obtained for the same 
condition. Once again, the primary features: splitting of the 
beam into two parts, strong scattering in the downstream 
direction and a weaker scattering in the upstream direction, 
are well represented. Similar to the shadowgraph situation, 
the experimental data also suffer from speckle noise which 
obliterates the weak fringes, predicted by the computational 
data. 

The primary shortcoming of the numerical calculation is 
the 2-D modeling of a three-dimensional (3-D) situation. The 
bow shock is axisymmetric; however, the present calculation 
assumes that the shock shape remains constant in the missing 
direction. A second source of error is in the assumption that 
the refractive index remains constant downstream of the 
shock. As the bow shock curves from its tip region, the fluid 
properties and, therefore, the refractive index changes. How- 
ever, the success of the numerical calculation in capturing all 


features of the light scattering phenomena indicate that the 
above concerns are of secondary importance. 

V. CONCLUSION 

The physical appearance of the light scattering effects in 
the two optical situations considered in this paper is very 
different. In the shadowgraph situation a wide central dark 
band and a few fringes are observed, while in a narrow laser 
beam and shock interaction situation long streaks in the 
farfield are observed. However, a closer look reveals many 
similarities. First, splitting of the main laser beam in the 
latter can be attributed to the shock shadow formation similar 
to that seen in the shadowgraph situation. Second, both situ- 
ations involve light propagating tangential to the shock sur- 
face. The bow shock, considered in the present work, is a 
3-D surface of revolution. In the shadowgraph image of Fig. 
2, light has passed through every part of this surface; how- 
ever, the shock shadow appears as a single line along the 
outer edge of the surface. Noticeably, only along this line, 
light rays are at a grazing incidence on the shock surface — 
identical to the situation when long light streaks are found to 
appear from a narrow laser beam. Third, upstream fringes 
and the upstream part of the light streaks are weaker than 
their downstream counterparts. 

Perhaps, a physical explanation unifying the two optical 
situations can be best provided from the Geometrical Theory 
of Diffraction. 16,17,9 According to this theory, electromag- 
netic diffraction of high frequency waves by straight edges 
gives rise to cylindrical diffraction edge waves. The total 
diffracted field, observed on the screen, is a sum of the inci- 
dent waves and the edge waves. In a shadowgraph situation 
the spatial extent of the incident field is very large and the 
edge waves interfere with this field to produce the fringe 
patterns seen upstream and downstream of the shock. On the 
other hand, the use of a laser beam of narrow spatial extent 
separates away the edge waves which then appear as the 
diverging light sheet from the point of its generation. A good 
example is the numerically calculated scattered light pattern 
shown in Fig. 12(b). The intensity distribution shows fringe 
formation only in the paraxial region where light from the 
undisturbed beam also appears. 

The effect of the refractive index jump, caused by a 
change in the free stream Mach number and total pressure, 
on the light scattering phenomena are studied in the present 
paper. Use of a different wavelength of light is also expected 
to change the scattering pattern, as in any other situation 
involving diffraction. A fourth parameter of interest is the 
thickness of the shock, which is estimated to be of the order 
of light wavelength 3 (light wavelength: 633 nm; shock thick- 
ness for the experimental range: 150-750 nm). Direct mea- 
surement of the shock thickness was not possible in the 
present experiment. Inclusion of this parameter in the nu- 
merical modeling requires complete solution of the Helm- 
holtz equation. Moreover, the experimental observation, that 
the light scattering phenomena are independent of the polar- 
ization of the incident field, indicates that the high-frequency 
approximation is valid for this situation. The shock, there- 
fore, is assumed to be of zero thickness for all analysis. 
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